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Abstract 

A unified model of dark energy and matter is presented using the modified variable 
Chaplygin gas for interacting dark energy in a non-flat universe. The two entities interact 
with each other non-gravitationally which involves a coupling constant. Due to dynamic 
interaction, the variation in this constant arises that henceforth changes the equations of 
state of these quantities. We have derived the effective equations of state corresponding to 
matter and dark energy in this interacting model. Moreover, the case of phantom energy is 
deduced by putting constraints on the parameters involved. 
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1 Introduction 



Hosts of results coming from the observations of WMAP [U [2] have convincingly shown the 
validity of standard Big Bang model with cosmological constant A, so called ACDM model. Sur- 
prisingly the energy density corresponding to A is two third of the critical density or ~ 0.7, 
the remaining part is due to baryonic and non-baryonic dark matter. The genesis of this energy, 
commonly called the dark energy, and its subsequent evolution poses serious problems to investi- 
gate. The dark energy is commonly represented by a barotropic equation of state (EoS) p = —p 
(or CO = —1). Numerous models have been proposed in the literature to explain the cosmic 
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accelerated expansion which include those based on homogeneous and time dependent scalar 
field termed as quintessence \3\, phantom energy having EoS a; < — 1 [1], Chaplygin gas (CG) 
with the EoS p = —X/p [5j, dissipative or viscous cosmology [6] and modified gravity models 
including Dvali-Gabadadze-Porrati (DGP) model [7J and higher order f{R) gravity theories [8]. 

Among all of these alternative theories, the Chaplygin gas efi^ectively explains the evolution of 
the universe right from the earlier matter dominated era (decelerating phase) to the late dark 
energy dominated era (accelerating phase). This CG was first introduced in the aerodynamical 
context [9] and has been largely investigated in recent years due to its adequate efficiency in 
interpolating the astrophysical data. 

The energy conservation equation p + 3H{p + p) = yields 



where X is a positive constant and Y is the constant of integration. Above, p and p are the 
pressure and energy densities of the CG. Thus small values of a(t) mimic dust phase while 
large values represent dark energy phase of the universe. The CG model favors a spatially flat 
universe with 95.4 % confldence level which agrees with the observational data of Sloan Digital 
Sky Survey (SDSS) and Supernova Legacy Survey (SNLS) [10]. It also possesses a property of 
giving accelerated expansion even if it gets coupled with other scalar fields like quintessence or 
dissipative matter fields [11]. Due to its this effectiveness, various generalizations of CG have 
been proposed in the literature (see e.g. [I2l[l3llllll5l[l6l[17lll8llia[20l[2l]). All these models 
based on the inhomogeneous Chaplygin gas offer unified picture of dark matter and dark energy 
(sometimes called dark matter-dark energy unification) [22]. But at the same time CG has a 
drawback as it produces oscillations or exponential blow up of dark matter power spectrum, 
which is inconsistent with the observations [23j . Similar results have been obtained in later 
generalizations of CG |24j . 

One of the long standing problems in cosmology for the past decade is the so-called "cosmic 
coincidence problem" which naively asks: If the universe evolves from the earlier quintessence 
{uj > —1) to the later phantom regime {to < —1), then why u = — 1 crossing is occurring at 
the present time. The astrophysical data seems to favor the evolving dark energy in which 
the EoS parameter ui changes at different cosmic epochs [25]. This problem can be rephrased 
alternatively as: In any cosmological model explaining dark energy, the ratio of matter to dark 
energy density is expected to decrease with the expansion of the universe. However the 
observations suggest that this ratio is of the order unity. It may implies the transfer of energy 
between these two entities to keep a delicate balance at the current time. The models in which 
this interaction is investigated are called the 'interacting dark energy models'. The interaction 
involves a coupling constant that determines the strength of this interaction. 

Earlier, it had been suggested that interacting models of CG with dark matter can effectively 
ameliorate the coincidence problem which is not possible in pure CG models \26\ 1271 [28]. More- 
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over, this problem has also been investigated in the context of dissipative CG and successfully 
explains the cosmic conundrum problem [11]. The interacting model effectively explains the 
phantom divide u = —1 i.e. the transition from a;>— ltow<— 1 [25^ \29\ I30j . It also rules 
out the possibility of the exotic phenomenon termed 'Big Rip' which is the infinite expansion 
of the universe in a finite time ^T]. It also yields stable stationary attractor solution of the 
Priedmann-Robertson- Walker (FRW) equations at late times. So if dark energy evolves then it 
would be worthwhile to find its effective EoS. This motivates us to investigate the more gener- 
alized form of interacting CG and to find the effective EoS of dark energy and matter in the 
respective model. The formalism adopted here is from Ref. |19j . 



2 Interacting modified variable Chaplygin gas 



We assume the background spacetime to be spatially homogeneous and isotropic represented by 
the FRW metric 

1 — kr^ 



ds' 



-dr + a\t) 



(2) 



Here a{t) is the scale factor and k = — 1,0,+1 represent spatially hyperbolic, flat and closed 
universes respectively. The corresponding Einstein field equation is 



3M2 



(pa + Pn 



(3) 



where = (SvrG) ^ is the modified Planck mass. Moreover, the energy conservation for our 
dynamical system is given by 

PA +/0m + 3-fr[(l +Wa)pa + Pm] = 0, (4) 

where we have taken the EoS for the dark energy pA = ^APA and for the matter part pm = 

(or LJm = 0). 

To derive the general form of the interaction between CG and matter, we start by considering 
the action of the scalar-tensor theory of gravitation given by 



'5T 



3f? 



1 



1 



62(0) 



(5) 



where is the Ricci scalar, (j) is the spatially homogeneous scalar field, represents the matter 
Lagrangian, ^ is the collective degrees of freedom of matter and is the coupling function 
between matter and the scalar field. Using the action Sst^ the interaction term Q can be written 
as [261 [32] 

Q = PmH— , (6) 

a ma 



3 



where 9' = e(a)(3'^™-^)/2^ Letting 



e'(a) = Goe' 



3 rh2(£l!l+£A)dlna 

J V Pm / 



Substitution of Eq. (7) in (6) yields 



Q = 3Hb'^{pA + p„ 



(7) 



(8) 



which is the general form of the interaction between CG and matter. Here b"^ is the dimen- 
sionless coupling constant (also called 'transfer strength') for the interaction. The positive 
is responsible for the transition from CG to matter and vice versa for negative b^. Sometimes 
this constant is taken in the range [0, 1] (see Ref. [26] )• Note that if 6^ = then it represents 
the non-interacting FRW model while b^ = 1 yields complete transfer of energy from CG to 
matter. Recently, it is reported that this interaction is observed in the Abell cluster A586 show- 
ing a transition of dark energy into dark matter and vice versa [33j . Therefore the theoretical 
interacting models are phenomenologically consistent with the observations. 

Let us assume an interaction Q = Tp\ between CG and dark matter, where T is decay rate. 
Then Eq. (4) can be divided into two parts corresponding to pA and pm as 



PA + 3H{1 + u;a)pa = -Q, 

Pm + 3Hpm = Q, 



(9) 
(10) 



respectively. Eqs. (9) and (10) show that the energy conservation for dark energy and matter 
would not hold independently if there is an interaction between them but hold for the total 
interacting system as manifested through Eq. (4). Further, we define the density ratio as 
r^n = To study how this density ratio evolves with time, we differentiate with respect to 
t and obtain 

Pm I Pm PA 



dt 

Using Eqs. (9) and (10) in (11), we get 



Pm 
PA 



UJA 



Pn 



PA 



3H 



Defining the effective EoS for dark energy and matter as [34] 



to 



eff 



i^A + 



3H' 



eff 



-1 r 

r„, 3H ' 



(11) 



(12) 



(13) 



which involve the contribution from the interaction between matter and dark energy. Using Eq. 
(13) in Eqs. (9) and (10), we get 



PA + 3Hil+ulf^)pA 
Pn, + 3H{1 + ujtl^)pm 



0, 
0. 



(14) 
(15) 
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From the standard FRW model, the density parameters corresponding to matter, dark energy 
and curvature are defined as 

" p,r~3Hm^' ^^^> 

^ _ PA _ PA (.rJ^ 

= (18) 

The parameter 0^ represents the contribution in the total energy density from the spatial cur- 
vature and is constrained as —0.0175 < Qk < 0-0085 with 95% confidence level by current ob- 
servations [35]. It is shown that a non-zero positive curvature parameter k allows for a bounce, 
thereby preventing the cosmic singularities without violating the null energy condition p+p >0 
[36]. 

The density parameters in Eqs. (16), (17) and (18) are related as 

+ = 1 + S^fe. (19) 
Using definition of and Eq. (19), we can write 

1 + — riyv 

In our generalized model, we choose the modified variable Chaplygin gas (MVG) given by 



(20) 



PA = ApA-^, B{a) = Boa-^. (21) 
Pa 

Here A, Bq and n are positive constants and < a < 1. Notice that for 5 = 0, Eq. (21) 
represents the barotropic EoS (or barotropic fluid). In general, the barotropic EoS p = Ap 
represents various types of mediums. For instance, A = —1 gives cosmological constant or de 
Sitter vacuum; for A = —2/3, it corresponds to domain walls; if ^ = —1/3 it gives cosmic 
strings; while ^ = corresponds to dust or matter; A = 1/3 yields the EoS for relativistic gas 
like photon radiation; A = 2/3 gives the perfect gas and A = 1 represents the ultra-stiff matter 
[37]. If B{a) = Bo in Eq. (21), it gives EoS of modified CG. Again if A = and a = 1, it 
represents the usual CG. Recently it is deduced using the latest supernovae data that models 
with a > 1 are also possible [38]. Notice that for ^ = 0, Eq. (21) yields a fluid with negative 
pressure. Such a fluid is generally characterized in the quintessence regime. 

The modified form of the CG is also phenomenologically motivated and can explain the flat 
rotational curves of galaxies [39]. The galactic rotational velocity is related with the MVG 
parameter AasVc = \plA and the density parameter p is related to the radial size of the galaxy 
as p = ^^^^^ . At high densities, the first term in MVG dominates and gives the flat rotational 
curve that is consistent with with observations for galaxies. The parameter A varies from galaxy 
to galaxy due to variations in 14 . 
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The density evolution of MVG is given by 

PA = (Aa-" + Ca-")TT^, (22) 

where s = 3(1 + a)(l + A), A = '^^^^"1^° and C is the constant of integration. Differentiating 
Eq. (22) we get 

PA = -— ^[(Aa"" + Ca-')^{Ana-'' + Csa"^)] (23) 
l + a 



Further, choosing the decay rate to be 

r = 36^77(1 + r„). (24) 
Substituting Eq. (24) in Eq. (9) and using Q = 3b^H{l + rm)pk, we obtain 

= ^J7T^[(^"'" + Ca-^)-i(Ana-" + Csa-')\ - b'l^^ _ i. (25) 
3(1 + a) i^A 



Using Eq. (20) and (24), we have 



r = ^b'H^^. (26) 
sza 



Thus using Eqs. (14), (25) and (26), we arrive at the effective EoS for dark energy as 

1 

3(1 + a) 



^ [(Aa-" + Ca"")-i(Ana-" + Csa-')] - 1. (27) 



Notice that the cosmological models with phantom energy arise when uj'^j^^ < —1 which is 
possible if the quantity in the square brackets in Eq. (27) is less then zero. Therefore, we have 
two possible cases: 

Case (1). If (Aa^" + Ca^^) > and {Ana~^ + Csa^^) < 0, then we have a constraint on a{t) 

z£l < a'-n < z£l^ (28) 
A An ^ ^ 

This is possible when ^ > and s > n. 
Case (2). If (Aa"" + Ca"") < and (Ana"" + Csa-") > 0, then we have 

(29) 



An A 
here a{t) will be restricted when > and s < n. 

Thus the phantom dark energy case is possible if scale factor a{t) is constrained as given in 
Eqs. (28) and (29). Also, Eq. (27) represents usual cosmological constant if the quantity in the 
square brackets vanishes identically. 

Also we note that the quantity in Eq. (22) must be positive so that a{t) > (-^) • Therefore, 
the minimum value of a{t) lies at 



A 



(30) 
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Therefore the universe can experience a bounce before reaching the singularity when a > amin 
which corresponds to the bouncing universe model. This bounce is possible due to our assump- 
tion of positive curvature (see Ref. [41] for a review on bouncing cosmologies). Moreover, the 
effective EoS of matter is determined from Eqs. (13), (20) and (24) as 

= -b'(l + ^). (31) 



Using Q,\ = 0.7 and = 0.3, we find from Eq. (31) 

oj^f = -3.336^. (32) 

Note that ojf^-^ depends only on the coupling constant. If there is no interaction then utf^-^ = 0. 
In particular if < 6^ < 1, it implies —3.33 < oo^^ — 0- Therefore in the interacting model, 
matter can also behave as a fluid with negative pressure which can cause the cosmic acceleration. 



3 Conclusion 



In summary, we have determined the effective equations of state for matter and dark energy in 
the generalized model of interacting Chaplygin gas. In particular for dark energy, the modified 
variable Chaplygin gas equation of state is used. Since dark energy interacts with matter then 
due to energy exchange, the EoS of the both the quantities will be modified. The effective EoS 
of Chaplygin gas oj'j^'^ also describes the EoS of phantom energy if the scale factor is constrained 
by either Eq. (28) or (29). As discussed, the interaction is observed in the near galactic cluster 
Abell A586, we may hope to determine the effective EoS of both matter and dark energy as 
well through future astrophysical observations that will finally put further constraints on the 
parameters in the model. Also it is shown that a universe filled with the Chaplygin gas can 
bounce from the singularity if a{t) is greater then amin, where the later describes the minimum 
value of scale factor at the bounce. Finally, it is shown that matter can also generate negative 
pressure along with dark energy to cause cosmic expansion. 
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